LAYERING IN THE ISING MODEL 

KENNETH S. ALEXANDER, FRANCOIS DUNLOP, 
AND SALVADOR MIRACLE-SOLE 

Abstract. We consider the three-dimensional Ising model in a half-space 
with a boundary field (no bulk field). We compute the low-temperature 
expansion of layering transition lines. 



1. Introduction and results 

We consider the Ising model in the half-space Z\ C Z 3 , with spins <7j = ±1, 
i G Z\ = {(11,12,13), 1-3 > 1}. The value —1 of the spin is associated with 
component or species A of a mixture and the value +1 is associated with 
component or species B, while the other half-space {i% < 0} represents a fixed 
given substrate or wall W , made of a third component or species. The formal 
Hamiltonian is 

(1.1) H ABW = J AB J2 - ^ + X> - + J WB £(1 + °i) 

<i,j> 43 = 1 13=1 

with energy contributions 2Jab, %Jwa, 2Jwb associated respectively to pairs 
of nearest neighbors AB, WA, WB. In the first sum, < i,j > are nearest 
neighbors in Z+. A wetting transition may occur when the bulk phase is B 
(or i?-rich) but the wall prefers A: J WA < J WB . 

At zero temperature, a macroscopic film of A will separate the wall from 
the bulk phase if Jwa + Jab < Jwb- One says that the wall is "completely 
wet" by phase A. Raising the temperature will favor the presence of a film, 
because the AB interface brings entropy. Therefore, at positive temperature, 
a film of A will always be present if Jwa + Jab < Jwb- There is no wetting 
transition, only complete wetting. 

On the other hand, if Jwa + Jab > Jwb, at zero temperature no A is 
present, and at low temperature the wall will be only partially wet by phase 
A. The density of B tends exponentially fast to the bulk density of B as a 
function of the distance to the wall. Raising the temperature now may produce 
a transition from partial to complete wetting: this is the wetting transition 
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predicted by Cahn [S] on the basis of critical exponents, and then confirmed 
by numerical and real experiments. 

The existence of the wetting transition has been proved mathematically in 
the two-dimensional Ising model pQ, but not in the three-dimensional Ising 
model. Let us simplify the notation to J = Jab and K = Jwb — Jwa, with 

(1.2) J > 0, < K < J. 

Let denote the +/— interface tension, defined for the Ising model in the 
full space Z 3 , without wall, with Hamiltonian equal to the first term of (11. ip . 
Frohlich and Pfister (see formula (2.20) and Fig. 2 in [8]) have proven, among 
other things: 

(1.3) K < ir ± Partial wetting. 

This is a non-perturbative result, valid for all temperatures < T < T c . 

We shall consider only low temperatures, and perturbative arguments (not 
fully mathematically rigorous), indicating that the partial wetting range is 
slightly wider than (11.31) . and includes first order layering transitions, as we 
now explain. Consider the model in a box A C Z\, with bottom layer at 
is = 1, and boundary condition a on the other five sides of the box. Let 
Ai = A R {is = 1}. The Hamiltonian (II. ip may be cast into the equivalent 
form 

(1.4) H K {a K \a) = -2J\A l \ + J £ (1 - a^) + Kj^i 1 + 

<i,j>nA^0 «3=1 

In the first sum, i,j are nearest neighbors in Z\ (so neither i nor j is in 
the wall), and <7j or aj should be replaced by <jj or Oj wherever % ^ A or 
j ^ A. In the second sum, i e A. The constant term in front is a convenient 
normalization. Boundary condition n, with n — 0,1,2, ... , is associated with 
the configuration n in Z\, given by 

(1.5) Oi = —1 if is < n, = +1 if is > n. 

A possible scenario for the wetting transition is as follows (see Fig. 1): Let 

< K < J with J — K small. At T = we have configuration 0, and for 
small T, we are close to configuration 0, call it state 0: in the thermodynamic 
limit, the probability that at a given i the spin a, differs from <jj, defined by 
(ll.5p with n = 0, is small. State n is defined similarly from configuration n, 
for any n. As the temperature is raised, a first order transition will occur, 
from state to state 1, then as the temperature is raised further, from state 

1 to state 2, and so on. The level of the stable state n goes to infinity as the 
temperature approaches the wetting transition temperature, which in this case 
is strictly below the roughening temperature. This scenario, with a sequence 
of first order layering transitions leading to the wetting transition, is part of 
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the general picture which emerged based upon various physical heuristics and 
Monte-carlo simulations (see [U [T2] and references therein.) 
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Fig. 1. Layering transition lines near T = 0. Dotted line 
shows a path from partial to complete wetting. 

Let t = e~ 4/3J < 1 and u = 2(3{J - K) = 0(t 2 ). Note that each factor 
of t corresponds to two plaquettes of the interface. We find the following 



approximation to the coexistence (first order transition) lines starting from 



(t = 0,u = 


0): 






0/1 


u = 


-ln(l-t 2 )+t 3 + C(t 4 ) 




1/2 


u = 


-ln(l -t 2 ) - t 3 + 5i 4 + 0(t 5 ) 




2/3 


u = 


- ln(l - t 2 ) - t 3 + At 4 - 4t 5 + 0(t 6 ) 


3/4 


u = 


-ln(l-t 2 )-t 3 + 4t 4 -6t 5 + ft 6 


+ 0{t r ) 


4/5 


u = 


- ln(l - t 2 ) - t 3 + At 4 - 6t 5 + ft 6 


- 51t 7 + 0{f) 


5/6 


u = 


- ln(l - t 2 ) -t 3 + At 4 - 6t 5 + ft 6 


- 53t 7 + 162t 8 + Oif 


6/7 


u = 


- ln(l - t 2 ) -t 3 + At 4 - 6t 5 + ft 6 


- 53t 7 + 160t 8 






+ (B 9 + 2)t 9 + 0(t 10 ) 




7/8 


u = 


-\n(l-t 2 )-t 3 + At 4 -6t 5 + ft 6 


- 53t 7 + 160t 8 


(1.6) 




+ B 9 t 9 + 0(t 10 ) 





Here B$ is a constant which we do not calculate, but we show it is the same for 
all interface heights n > 6. The analogous statement applies to the calculated 
coefficients as well, for example, the coefficient of t 4 is 4 for all n > 2. This is 
a result of the cancellation of all terms proportional to n, n 2 , etc. in the low 
temperature expansion of the increment of surface free energy from nton + 1, 
up to the given orders in t. We are unable to determine a systematic way in 
which this cancellation occurs, but we anticipate its validity for all orders in 
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t. The consequence is that each successive transition line requires one more 
order in t to discern it. 

The phases 0, 1, 2, 3, 4, 5, 6, 7 are predicted to be stable between the 
respective transition lines. In particular phase should be stable for u > t 2 + 
t 3 + 0(t 4 ). For comparison, (11.3p gives partial wetting for u > 2t 2 + 4t 3 + 0(t 4 ). 
Basuev [3] has given such equations for coexistence of the phases 0,1,2 with 
1,2,3 respectively. 

Naturally, more is known in the SOS approximation, and in that context 
full mathematical rigor is possible, see [5J [2]. The low-temperature expansions 
of the Ising model and the corresponding SOS model agree only up to and 
including order t 2 , which is of little help for (11.60 . Order t 3 corresponds to a 
domino excitation of the interface, same in Ising and SOS, but also to a unit 
cube bubble, present only in the Ising model. 

The stability range of phase n appears to be of width approximately 2t n+2 
in the variable u. This is the same for Ising and SOS, and is the result of a 
double leg interface excitation reaching the wall (see Fig. 7). 

The n/n + 1 coexistence lines are expected to converge as n — > oo to a part 
of the wetting transition line. Therefore the low-temperature expansion of the 
n/n+1 coexistence lines for all n would give the low-temperature expansion 
of the wetting transition line. 

The derivation of the 2/3, 3/4, 4/5, 5/6, 6/7 transition lines is given in 
Section [21 except for the recursion diagrams, which are displayed and explained 
in Section HI The special features of the 0/1 and 1/2 transition lines are given 
in Section [3j Diagrams for the 7/8 transition line are postponed to Section [5j 



Let us consider a finite volume and boundary condition n, with n > 1 
for definiteness. The ground state is (11.51) . with a flat interface at height 
n + |, denoted /„. At positive temperature, bubbles and interface excitations 
will appear. If state n is stable, or if the statistical ensemble is restricted 
by a condition forbidding large fluctuations, the gas of bubbles and interface 
excitations should be diluted, and the corresponding dilute gas expansion is 
expected to give exact asymptotics for low temperatures. The corresponding 
partition function is 



2. LOW TEMPERATURE EXPANSION 



(2.1) 




where (3 = 1/kT is the inverse temperature and the ' indicates that summation 
is over a restricted ensemble corresponding to state n. The associated surface 
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free energy density (times ft) will be denoted /„, so that 



A 



(2-2) f n -fn+i= lim -J-flog-S- 

We are going to compute the leading terms up to some order for /„ — f n +i, so 
as to obtain ( 11.6f) . 

Bubbles and interface excitations will be called contours, or also polymers, 
and will be denoted 7. They are defined as boundaries of maximal connected 
sets of points where the spin differs from its ground state value in the corre- 
sponding restricted ensemble. A set of points is connected if any two points 
can be connected by a path of nearest neighbor bonds in the set. The bound- 
ary of a set of points is a set of plaquettes. A contour need not be connected. 
Interface excitations are distinguished by the property of sharing at least one 
plaquette with I n . A bubble crossing I n without sharing a plaquette is not an 
interface excitation. 

The low-temperture polymer expansion starts with 

(2.3) ^ = e B|Al " (B) EII^) 

{7} 7 

where {7} is a compatible family of contours, and ^(7) is the weight of a 
contour, 

(2.4) ^(7) = ^H-l7ni«l e u|7n{*=i}| 

where | ■ | is the number of plaquettes in 7 or in jC\I n or in 7IH {z = -}. A family 
is compatible if any pair of contours in the family is compatible. Two contours 
are compatible if their interiors are disjoint and they share no plaquette. In 
view of (12. 4p . we will represent an interface excitation with plaquettes in I n 
removed (see Fig. 2-7 below), but when deciding compatibility, it must be 
remembered that these plaquettes do belong to the interface excitation. 

As the interaction between contours is a two-body interaction — compati- 
bility is decided two by two — the general theory of polymer expansion (see 
e.g. EH HQ ) gives, from (Q, 



(2.5) M^ A ) = ;>> 



where u is a cluster or family of contours, with contour 7 repeated n 7 times, 
and 

(2.6) ^ T M=n(^T^) n7 )E(- i ) / 

7£w 7 ' G 
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where the sum over G is over connected graphs on the cluster, and I is the 
number of edges in G. An edge may exist between 7 and 7' if and only if 7 
and 7' are incompatible. 

For the expansion of r 1 * 1 , interface excitations were expanded in terms of 
walls and ceilings by Dobrushin [7j, who proved convergence of the resulting 
expansion. For the SOS approximation of the present wetting model, a two- 
scale convergent expansion was used in [2]. Here we consider only the finite 
volume expansion and the formal infinite volume series, which is why our 
derivation of (jl.fip is not fully rigorous. 

All the clusters in (12.51) lie within A. For a cluster which contains an interface 
excitation, we write uj G I n . For a cluster of bubbles only, compatible (i.e. not 
sharing a plaquette) with I n , we write to ~ /„. For a cluster which reaches the 
bottom {i% = 1/2}, we write uj G W, otherwise uj ~ W. We write Wn for the 
top boundary {23 = N + |} of A. All clusters uj C A are compatible with the 
top boundary; we write u ~ Wn- Then 



where 



(2.8) ^0(7) = t^- hnInl , ^(7) = t5lTl e «l7n{*=i}| ; ^( 7 ) = £ |H 



The first term in (12.71) depends explicitly upon n. The sums consist of clus- 
ters uj C A, but in order to extract the n-dependent part of the following 
three terms, it is convenient to relax this condition into uj H A 7^ 0, allow- 
ing "boundary-overlapping" clusters which overlap W or Wn- In this context 
the notations uj ~ W, uj G W and W ~ Wn apply only to clusters which do 
not overlap W and Wn respectively. Then applying inclusion-exclusion to the 
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summation conditions, the last three sums in (12.71) become 

E ¥ T (u) = E M - E ( w ) - E H + E M 

wG^ni UJGin u£ln, uj^I n , uj^I n , 

oj&W,W n uj^W uj^W n lu^W,W n 

E = E ( w ) - E m - E ^) + E 

oEIC, wGVF w6W, w6W, oEW, 

E ^ T( ^) = E ^ (w) - E ^ (w) - E ^ n - E ^2 m 

(2.9) E ¥#M+ E ^M- E ^H- 

Note that the sums from (12 .7p . on the left side in (12.91) . are not affected by 
the relaxation from ucAtowflA^I. Terms with uo ^ I n ,ui ^ Wn or 
uj tfi W, Wn or u> ^ W, Wn, ^ I n are negligible in the thermodynamic limit 
and will be omitted in the sequel. This is the meaning of ~ instead of = below. 
Apart from these negligible terms, only one sum on the right side each of the 
three equalities in (12. 9p actually depends upon n. Therefore 
(2.10) 

log(^n) - E ^ T M" E VoH - E E ^H+indep. of n. 

L0£l n ,W w€W, u^I n , 

uj^W ujtfiW 

In order to compare and Z^ +1 using translation invariance, the wall W will 
be denoted Wo, and W-\ will denote a wall translated vertically by —1. The 
following is immediate from (I2.10p . 

Proposition 1: For n > 1, in the limit of a box A of height N — > oo, 



iog(^ A /^ +1 )= E E ^M- E 

UEln,W uEl n+1 ,W 

-( e <£{<»)- E ^m) + (E^h- E 



1 UJ 



jew, uiew, ui^w, um^w, 

u^I n U^I n+ i u^In u "i iI n + l 

<"~Vfl u>~I„ 

(2.11) 

In terms of surface free energy densities, anticipating a leading term t 2n , this 
can be written as 

(2.12) r 2n (f n+1 - /„) = A n (u) - A n (0) - t 2 A n+1 (u) - B n (u) + B™(0) 

where each of the five terms is defined by the corresponding term in (12. lip . 
We can simplify B™(0) as follows. The terms in B^(0) correspond to clusters 
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of bubbles only, and the set of such clusters may be divided into equiva- 
lence classes consisting of clusters which are vertical translates of one another. 
Within each equivalence class there is a unique special bubble u satisfying 
uj G W . For a given equivalence class, the number of terms from that class in 
the first sum in 5^°(0) is the number of heights k > n for which the special 
bubble has a horizontal plaquette at height k + ~, and similarly for the second 
sum, but with heights k > n + 1. Hence the net number of terms in B^(0) 
from the equivalence class, counted with +/— sign, is 1 if the special bubble u 
has a horizontal plaquette at height n+| (that is, if u 7^ I n ), and otherwise. 
It follows that 

(2.13) t 2n B™(0) = 
Since 

t 2n B n ( U ) = («) - E 

uew, ci>ew, 

and since cpf = <p\ for u = 0, we have B~(0) = B n (0) + ^^(O) so that 

t 2U (fn+l — fn) 

(2.14) = A n (u) - A n (0) - t 2 A n+1 (u) - (B n (u) - B n (Q) - t 2 B™ +1 (0)). 
The u dependence may be written as 

(2.15) A n {u) = e u P n + e 2u Q n + e 3u R n + e iu S n + e 5u T n + e Qu U n + ... 

(2.16) B n (u) = e u P n + e 2u Q n + e 3u R n + e 4u S n + e 5u f n + e 6u U n + ... 

For n > 3 we have P n = 0(1) corresponding to interface fluctuations placing 
a single plaquette on the wall, and similarly Q n = 0(t 2 ), R n = 0(t 4 ), S n = 
0(t 5 ), T n = 0(t 7 ), U n = 0(t s ). Relative to these, P n , Q n ,R n ,S n have an 
extra factor t at leading order. The remainder in f)2.15p . (12.161) is 0(t 10 ). For 
n = 2 we have P 2 = 0(1), Q 2 = 0(t 2 ), R 2 = 0(t 4 ), S 2 = 0(t 4 ), T 2 = 0(t 6 ), 
U 2 = 0(t 6 ), while P 2 ,Q 2 , R 2 , S 2 are of the same order as for n > 3. The 
remainder in fLTTBD for n = 2 is £>(t 8 ), but in fLTTBD it is still 0(t 10 ). 

Let Q n = Qn + Q n and R n = R n + R n + R n , where the upper index 1, 2, 3 is 
the number of polymers (in the cluster) touching the wall, so that Q n = 0(t 2 ), 
Q 2 n = 0(t 3 ), etc. We are going to expand (12 . 1 2 j) up to order t 9 , requiring A n +i 
up to order t 7 , using recursion in n. 



00 

= t 2nJ2t 2m B n+m (0). 
m=0 
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(2.17) 

Pn+i = P n + 2Qn + 3P n + AS n + 5T n + 6U n -t(P n + 2Q n + 3R n + AS n ) 
+ 0(t 7 ) 

Qn+i = (4t 2 - 4t 3 )P„ + (t + 6t 2 - 7t 3 )Q* + (8t 2 - 8t 3 )Q 2 

+ 2tR 1 n + 9t 2 P n + tR 2 n + AtS n + 0(t 7 ) 
Ql+i = ("5t 3 + 5t 4 )P n + (-10t 3 + 10t 4 )g^ + t 2 Q 2 n - l2t 3 Q 2 n + 0{t 7 ) 
Q n+1 = (At 2 - 9t 3 + 5t 4 )P n + {t + 6t 2 - 17t 3 + lOt^Ql + (9t 2 - 20t 3 )Q 2 n 

+ 2tR\ + 9t 2 R n + tR 2 n + AtS n + 0(t 7 ) 
R} n+1 = (18t 4 - 18t 5 )P n + (6t 3 + 2U 4 )Ql + t 2 R n + U 2 S n + 0{t 7 ) 
R 2 n+1 = (-48i 5 + 48t 6 )P n - 8t 4 Q^ + 0(t 7 ) 
R 3 n+l = 31t 6 P n + Oif) 

R n+1 = (18t 4 - 66t 5 + 79t 6 )P n + (6t 3 + 16t 4 )Q^ + t 2 R n + At 2 S n + 0{t 7 ) 
S n+1 = (At 5 + 60t 6 )P n + 2t 4 Q n + 0(t 7 ), 



and the same recursion relations for P n , Q n , R n , S n , with an error 0(t 8 ). The 
recursion relations (12. 17j) have been found with the help of diagrams, see Sec- 
tion n 

Solving these recursion relations for formal power series in t requires as 
input P n , or P n , for all n, to the required order. Indeed the order obtained 
in the output P n+ i or P n+ i is the same as in the input P n or P n , so that 
the recursion formula does not help. On the other hand, if the power series 
expansion for P n or P n is obtained by other methods, up to the required order, 
for all n, then the initial condition, at n = 2, given by Q\ = At 2 + 2t 2 + 0(t 3 ), 
Ql = -5t 3 + 0{t 4 ), R\ = 0{t 4 ), R 2 = 0{t 5 ), S 2 = 0{t A ), or Q\ = At 3 + 0{t A ), 
Ql = -5t 4 + 0{t 5 ), R\ = 18t 5 + 0{t 6 ), together with P n or P n , will give one 
more order in t with each recursion step. The recursion equation giving P n +i 
or P n +i may be checked at the end for consistency. The final result of this for 
P n , Q n , R n , S n or P n , Q n , R n , S n is given as "first excitations": 
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First excitations in A n (u): 
(2.18) 

P n = 1 - (n - 5)t + c n t 2 - a n t 3 + d n t A + 0(t 5 ) , n > 5 

Q l n = A[t 2 - (n - 6)t 3 + (c n _i + l)t A - (a„_i + c„_i - c n _ 2 + 6n - 41)t 5 

+ (rf n _i + a n _i - a n _ 2 + 5c n _ 2 + c n _ 3 + 2n + C)t 6 ] + 2t n + 0(t 7 ) ,n>6 
Ql = -5t 3 + 5(n - 5)t 4 - (5c n + C)t 5 + C(t 6 ) , n > 3 
Rl = 18t 4 - (18n - 114)t 5 + (18c n - 24n + C)t 6 + 0(t 7 ) , n > 4 

i?2 = _ 48t 5 + ( 4gn + ^^6 + ^7j ; n > 3 

5 n = 4t 5 - (4n + C)t 6 + C(t 7 ) , n > 3 
T n = Ct 7 + Oif ) , n > 3, 

with 

(2.19) c n +4(n-2) + 16, 

(2.20) a n = ^ + 12 2 ^ " 10n _ 48 ' 

(2.21) d n = ~ ^ + 20^ 3 X ) + ^ X ) + 54 ™ + a 

In (12.181) and in what follows, C is a generic constant, not depending on n 
and different at different appearances, which we do not calculate or use. The 
expansion for P n is valid for n = 3 with two orders less (that is, 0(t 3 ) instead 
of 0(t 5 )), and for n = 2 with three orders less, and for n = 1 with four 
orders less. It is obtained by listing diagrams — see below. The results for 
Q\, . . . S n in (12.181) follow from the result for P n using (I2.17p . One can start 
the induction from n = 1 with (12.171) adjusted for n = 1, or from n = 2 
with Ql = At 2 + 2t 2 + 0(t 3 ), Q 2 = -5t 3 + 0(t 4 ), R\ = 0(t 4 ), R 2 = 0(t 5 ), 
5*2 = 0(t 4 ). The expansion for Q\ is valid for n = 4 with one order less, and 
for n = 3 with two orders less, and for n = 2 with three orders less. The 
expansion for R^ is valid for n = 3 with one order less. 

The result for P n is displayed in Figs 2-6. Formula (12.191) for c n was obtained 
using Fig. 2 and Fig. 3 with (12. 6p . The factor 6 for the last diagram in Fig. 
3 is: one incompatible unit cube upward interface excitation, as drawn, and 
five incompatible unit cube downward interface excitations. Formula (12.201) 
for a n was obtained using Fig. 4, 5, 6 with (I2.6p . The factor 5(n — 2) for the 
one before last diagram in Fig. 4 is: one incompatible unit cube bubble in the 
interface leg at height 2, . . . , n — 1, as drawn, and four incompatible unit cube 
bubbles adjacent to the leg at height 2, . . . ,n — 1, and similarly for the last 
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diagram in Fig. 4. Formula ( 12.211) was obtained using the diagrams in Fig. 13 
and Fig. 14. 



-At -(n-l)i +( n 2 1 )t 2 -4(n-l)< 2 +4(n-2)t 2 +4(n - 2)t 2 
Fig. 2. P n : up to the t 2 -terms dependent on n. 



r 



-I2t 2 



-U 2 



+At 2 



-At 2 



-6t 2 



Fig. 3. P n : i 2 -terms independent of '. 



'n-i\ f 3 



t A +4 



)i 3 



(Y)* 3 -5(n - 2)t 3 -5(7i - 2)i 3 



Fig. 4. P n ; t 3 -terms, cubic, quadratic (all) or linear (continued on next two Figs.) in n. 



r 



-12(re- l)t 3 



-6(n-l)t 3 -A(n-2)t 3 -4(n - l)i 3 +6(n - 1)< 3 
Fig. 5. P n : i 3 -terms, analog of i 2 terms on Fig. 3. 



12 



K.S. ALEXANDER, F. DUNLOP, AND S. MIRACLE-SOLE 



+(16(n - 2) - 4)t 3 +4(n - 3)i 3 +(16(n - 2) - 4)t 3 +4(n - 3)t 3 

Fig. 6. P n : i 3 -terms, linear in n (continued from previous Figs.). 



The leading terms up to t 3 and the double leg in Q n = Q\ + Q\ are shown on 
Fig. 7. 




At 2 -4(n - l)t 3 +16i 3 +4t 3 -5i 3 

Fig. 7. Q™ = + Qn : U P to order i 3 , and double leg. 



+2t Tl 



Formulas RHM and (1X20]) for c n and consistent with the recursion 

relations, notably the equation giving P n+ i not used so far, implying 



(2.22) 



c n+ i = c n + n + 3 

a n+ i = a n + c n + 8n — 30 . 



For later purposes we note that 



(2.23) a n - a n -i + 2c n - 4c„_i + c n _ 2 = 9n - 35 



l)(n - 2) + lln - 32 
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Putting together ( 12. 15ft . ( 12.171) and assuming u = 0{t 2 ) gives for n > 3 
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(2.24) 

A n (u)-A n (0)-t 2 A n+1 (u) 

= (e u - 1 - e u t 2 + e u t 3 - 4e 2 V + 9e 2u t 5 

- 5e 2u t 6 - 18e 3u t 6 + 66e 3u t 7 - Ae 4u t 7 - 139t 8 + Ct 9 )P n 
+ (e 2u - 1 - 2e u t 2 + 2e u t 3 )Q n + (-e 2u {t 3 + 6t 4 ) + lit 5 - 28£ 6 + Ct 7 )Q\ 
+ (-9t 4 + 20t 5 + Ct 6 )Q 2 n + (e 3u - 1 - 3e u t 2 + 3e u t 3 )R n 
- (2t 3 + 10t 4 + Ct 5 )Rl - {t 3 + Ct 4 )R 2 n 
+ (e 4u - 1 - Ae u t 2 + 4e u t 3 + Ct 4 )S n - {At 3 + Ct 4 )S n 
+ (e 5u - l-5e u t 2 )T n + 0{t 10 ). 

Contributions from U n have been absorbed into 0(t 10 ), thanks to u — 0(t 2 ). 
For n = 2, (jg]2lj) is valid up to order t 6 , with an error 0{t 7 ). By (|2TT8|) . 
in each of the expansions P n , Q\, Q 2 n , etc., n-dependent terms only appear at 
one or more orders less in t than the largest-order term, and when (12.241) is 
multiplied out, the unspecified constants C only appear at order t 9 or less. 
Therefore the constants C appear in n-dependent terms only at order t 10 or 
less, so that, while the constants C are relevant to the value of B$ in ( II. 6p . 
they are not relevant to establishing that Bg is n- independent. The constants 
C thus play the role of placeholders, permitting the 0(t 10 ) error term which 
allows analysis of the dependence of Bg on n. 



First excitations in B n (u): 

P n = t-(n- l)t 2 + c n t 3 - ~a n t A + d n t 5 + Oif) , n > 4 

Qn = 4[t 3 - (n - 2)t 4 + (cn_i + 7)t 5 + (-a n _! - c n _i + c„_ 2 - 6n + 17)t 6 ] 

+ 2t n+2 + 0(t 7 ) , n > 4 
Q2 = _ 5t 4 + 5 ( n _ ^ t 5 + ^6^ ^ n > 2 

= I8t 5 - (18n - 42)t 6 + Oit 7 ) , n > 2 
Rl = -48t 6 + C(t 7 ) , n > 2, 
(2.25) 
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with 



(2.26) c n 

(2.27) a n 

(2.28) d n 



The expansion for P n is valid for n = 3 with one order less, and for n = 2 
with two orders less, and for n = 1 with three orders less. The expansion 
for Q n is valid for n = 3 with one order less, and for n = 2 with two orders 
less. Formula (12.261) was obtained using Fig. 8. The last two diagrams in 
Fig. 8 belong to the second term inside the parentheses, 4th term in (12. lip , 
defining P n . Formula (I2.27P was obtained using the last diagram in Fig. 8 
and all diagrams in Fig. 9. Formula (12.281) was obtained using the diagrams 
in Fig. 15. 



n — 1 
2 

n — 1 
3 

n — 1 
4 



+ 8n - 13 



16 
24 



n — 1 
2 

n — 1 
3 



+ 5n + 2 



79 



n — 1 
2 



- 31n + 75. 



< -(n- l)i 2 + (' l 2 1 )i 3 +4(n-2)i 3 +4(n - l)i 3 -< 3 +(n - l)t 4 
Fig. 8. P n : i 3 -terms, and t -term incompatible with / n +i. 



r - i 



- ( n - 1 ) t 4 -4(n - l) 2 i 4 -8 ( n 2 2 ) t 4 -5(n - l)t 4 -5(n - 2)i 4 +4(n - 3)< 4 +4(n - 2)< 4 
Fig. 9. P n : t 4 -terms, other than Fig. 8. 
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These c n and a n are consistent with the recursion relations, which imply 

Cn+i = c n + n + 7, 

(2.29) a n+ i = a n + c n + 8n + 2. 

Then 

c n - c n _i = n + 6, 

(2.30) 5 n+2 - a n = (n - 1) (n - 2) + 33n - 7, 

and from this we obtain 

(2.31) a n+2 - a n + 2c n - 4c n _i = 21n + 43. 

Putting together (12.161) and the analog of (12.171) for P n , Q n , etc., and assuming 
u = 0(t 2 ), gives for n > 2 

B n (u)-B n (0)-t 2 B n+1 (0) 

= (e u - 1 - t 2 + t 3 - 4t 4 + 9t 5 - 23t 6 + 62t 7 - 139t 8 )P n 
+ (e 2u -l-2t 2 + 2t 3 )Q n + (-t 3 - Qt 4 + lit 5 - 28t 6 )Q* 
- (9t 4 - 20t 5 )Q 2 n + (e 3u - 1 - 3t 2 + 3t 3 )i? n - (2t 3 + 10t 4 )i^ 

(2.32) -t 3 R 2 n + 0{t 10 ), 

while from (I2TT6D and (IZ25I) . 

t 4 5„ +2 (0) = t 5 - (n + l)t 6 + (5 n+2 + 4)t 7 - (& n+2 + 4n + 5)t 8 

(2.33) + (d n+2 + 45 n+1 + 5n + 4l)t 9 + 0(t 10 ), 

(2.34) t e B n+3 (0) = t 7 -(n + 2)t 8 + (d n+3 + 4)t 9 + 0(t 10 ), 

(2.35) t 8 B n+A (0) =t 9 + 0(t 10 ). 

Assumption: 

u = - ln(l - t 2 ) + 6 3 t 3 + ■ ■ ■ + + M 9 + 0(t 10 ) 

=t 2 + 6 3 t 3 + (6 4 + |)t 4 + M 5 + (6e + |)^ 6 + M 7 + (6 8 + \)t 8 + b 9 t 9 + 0(t 10 ) 
(2.36) 



Then, with 6 3 = 



— 1 and 64 = 4 where 63 and 64 don't appear explicitly, 
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e «(! _ t 2 ) - 1 + e"t 3 = (6 3 + l)t 3 + M 4 + (65 + l)t 5 + 0(f) 

= At 4 + (65 + l)t 5 + (h - |)t 6 + (67 + I)* 7 
+ (6 8 + 7)t 8 + (69 + C)t 9 + 0(t 10 ) 
e 2u -l- 2e u {t 2 - t 3 ) = 2{b 3 + l)t 3 + (26 4 + l)t 4 + C(t 5 ) 

= 9t 4 + 26 5 t 5 + (2b 6 + 10)t 6 + C(t 7 ) 
e 3u - 1 - 3e u (i 2 - t 3 ) = 3(63 + l)t 3 + 3(64 + l)t 4 + 0{t 5 ) 
(2.37) e 4u - 1 - 4e"(t 2 - t 3 ) = 4(6 3 + l)t 3 + (46 4 + 18)t 4 + 0{t 5 ) 



and 



e" - 1 - t 2 + t 3 = (& 3 + l)t 3 + (64 + l)r + C(t 5 



= 5t 4 + (65 - l)t b + (b 6 + ^)t b + (67 + 65 - 5)t 7 
+ (-&5 + &6 + &s + f)t 8 + 0(t 9 ) 
e 2 " - 1 - 2t 2 + 2t 3 = 2(63 + l)t 3 + (26 4 + 3)t 4 + Oif) 

= lit 4 + (26 5 - 4)t 5 + (26 6 + 22)t 6 + Oit 1 ) 

(2.38) e 3u - 1 - 3t 2 + 3t 3 = 18t 4 + Oif ). 
Then for n > 2, from fl2T2lj) . (ESTJ), (12351) . 

A» - A n (0) - t 2 A„ +1 (n) = [(63 + l)t 3 + (64 - 4)t 4 ]P„ + 0(t 5 ) 

(2.39) =(6 3 + l)t 3 + [64 - 4 - (63 + l)(n - 5)]t 4 + 0(t 5 ), 
while from fT2T32l) . fl2T33|) . fl2T38l) . fl2T25|) . 



B n {u) - t 2m B n+m (0) = [(63 + 1)^ 3 + (64 - 3)t 4 ] P n - + ^(t 6 



m=0 



(2.40) 



giving 



(63 + l)t 4 + [64 - 4 - (63 + l)(n - l)]t 5 + 0(t 6 ), 



(2.41) £~ 2ri (/ n+1 - f n ) = (63 + l)t 3 + [64 - 4 - (63 + l)(n - 4)]t 4 + 0(t 5 ) 
If 63 = -1, then from fl2T24j) . (12371) . fl2TT8|) . still for n > 2, 

A n (w) - A n (0) -t 2 A n+1 (w) 

= [(64 - 4)t 4 + (65 + 10)t 5 ] P n - t 3 Q n + 0(t 6 ) 

(2.42) = (64 - 4)t 4 + [65 + 6 - (64 - 4)(n - 5)]t 5 - 2t"+ 3 + C(t 6 ), 
giving 

(2.43) r 2n (/„ + i-/„) = (64-4)t 4 +[6 5 + 6-(64-4)(n-4)]t 5 -2r +3 + 0(t 6 
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A n (u) - A n (0) -t 2 A n+l (u) 

= [(6 6 + 10)t 5 + (6e - I - 31)£ 6 ]P n + 9t A Q n - (t 3 + Qt A )Q l n + 0(t 7 ) 

(2.44) = (& 5 + 6)t 5 + [& 6 - f - (& 5 + 6) (n - 5)] t 6 - 2t n+3 + 0{t 7 ) 

while from (12321) . (1238]) . (T2~25D . 

00 

B n (it) - ^ t 2m P„.+ m (0) = [t 4 + (65 + 8)t 5 ]P n - t 3 Q n -t 5 + (n + l)t 6 + 0(t 7 ) 

(2.45) = (65 + 6)t 6 + 0(t 7 ) 

giving 

(2.46) t- 2n (f n+1 -f n ) = (b 5 + 6)t 5 +[b 6 -f-(b 5 + 6)(n-A)]t e -2t n+3 + 0(t 7 ) 

If 6 5 = -6, then from (12T24D . (12371) . (12351) . now for n > 4, 

A n (w) - A n (0) -t 2 A n+1 (u) 

= [4t 5 + (6e - \ - 31)t 6 + (67 + 89)t 7 ]P n + [9t 4 - 12t 5 ]Q n 

- (t 3 + 6t 4 - 9t 5 )gi - 9t 4 Q2 _ 2t 3 Rl + G(t 8 ) 

= (b 6 - f )t 6 + [67 + 4(c n - c n _! - 3n) - (6e - f )(n - 5) + 85] t 7 

- 2t n+3 + 0(t 8 ) 

= (&e - f + [67 + 53 - (6e - f )(n - 5)]t 7 - 2t n+3 + 0(t 8 ), 

(2.47) 

while from (12321 . (12381) . (LT251) . 

S n (u) - B n (0) - t 2 B n+1 (0) 

= [t A + 2t 5 + (be - f )t e + (67 + 51)t 7 ]P„ + (lit 4 - lQt 5 )Q n 

+ (-t 3 - Qt A + llt 5 )Q l n - 9t A Q 2 n - 2t 3 R n + 0(f) 
= t 5 -(n + l)f + [c n + 2n + 6 6 - \]t 7 
(2.48) + [b 7 - ~a n + 2c n - 4c n _! - (6 6 + |)n + b 6 -f]t 8 + 0(t 9 ) 



18 K.S. ALEXANDER, F. DUNLOP, AND S. MIRACLE-SOLE 

so that, with (ET331 . fLTM|) . 



OO 

= (6e - f + [67 + ^6 + S n+2 - a n + 25 n - 45 n _! - (6 6 - |)n - f ]t 8 
+ 0(t 9 ) 

= (6e - f)t 7 + [67 + 53 - (&e - f )(n - l)]t 8 + (9(t 9 ), 

(2.49) 



giving 

(2.50) r 2n (f n+1 -f n ) = (6 6 -f )t 6 +(6 7 +53-(6 6 -f )(n-4))t 7 -2r +3 +G(t 8 ). 

If 6 6 = f , then from (12331) . fl237j) . fl2TT8|) . now for n > 5, 

- A n (0) -t 2 A n+1 (u) 
= [At 5 - 8t 6 + (67 + 89)t 7 + (& 8 - 258)t 8 ] P n + (-t 3 + 3t 4 - 3t 5 + 19t 6 )Q, 

+ 8t 5 Q* + (-2t 3 + 5t A )R l n - t 3 R 2 n - At 3 S n + 0(f) 
= (67 + 53)t 7 

+ [& 8 - 4(a n - a n _i + 2c n - 4c n __i + c„_ 2 ) 

- (& 7 + 53) (n - 5) + 36n - 300]t 8 - 2t n+3 + 0(t 9 ) 
= (67 + 53)t 7 + [6 8 - 160 - (67 + 53) (n - 5)]t 8 - 2r+ 3 + 0{t 9 ), 
(2.51) 



while (12.491) becomes 

00 

(2.52) B n (u) - t 2m B n+m (U) = (67 + 53)t 8 + G(t 9 ) 

m=0 

giving 
(2.53) 

t- 2 "(/n+i - /n) = (67 + 53)t 7 + [6 8 - 160 - (& 7 + 53) (ti - 4)]t 8 - 2t" +3 + 0{t 9 ). 
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Finally, if 67 = —53 then 
A n (u) - A n {0) -t 2 A n+l {u) 

= [4t 5 - 8t 6 + 36t 7 + (6 8 - 258)t 8 + (69 + C)t 9 ] P n 

+ (-t 3 + 3t 4 - 3t 5 + 19t 6 + Ct 7 )Ql + (8i 5 + Ct 6 )Q 2 n 
+ (-2t 3 + 5t 4 + Ct b )R l n - (t 3 + Ct^Rl - (At 3 + Ct 4 )S n + Oif) 
= (b 8 - 160)t 8 + [69 + 4{d n - d n -i) + 8a„ - 12a n _i - 4(a n _ x - a„_ 2 ) 

- 24c„_! - 8c n _ 2 - 4c n _ 3 - 92n + C]t 9 - 2t n+3 + 0(t 10 ) 
= {b 8 - 160)t 8 + [69 - C - {b s - 160) {n - 5)]t 9 - 2t n+3 + 0{t 10 ), 

(2.54) 

while from (12321 — (12351) . (I23B1 . (12321) . 

00 

S n (u)- ^ ^5,^(0) 

m=0 

= [t 4 + 2t 5 + 6t 6 - 2t 7 + (6 8 - 96)t 8 ] P n + (lit 4 - 16t 5 + 69t 6 )Q„ 

+ {-t 3 - 6t 4 + lit 5 - 28t 6 )gi - (9t 4 - 20t 5 )Q 2 n + {-2t 3 + 8t 4 )P 4 

- t 3 R 2 n + 0{t w ) 

= [bs + dn~ dn+2 ~ 2a„ + 4a n _i - c n+3 - 4c n+1 + 6c n + 20c„,_i - 4c n _ 2 

+ 87n + 130]t 9 + O(t 10 ) 
= (6 8 -26)t 9 + 0(t 10 ), 
(2.55) 

giving 
(2.56) 

r 2n {f n+1 - / n ) = (6 8 - 160)t 8 +[b 9 -C- {b 8 - 160) (n - 4)] t 8 - 2t n+3 + 0(t 9 ). 
Now, collecting (jCTj) . ([233]) . (12361) . (1230|) . (1233|) . (12361) gives: 

Proposition 2: The following are valid for sufficiently small t. 

• If 63 > -1, or 63 = -1, 64 > 4, then 

(2.57) /«+i-/n>0, n>2, 

and phases 3,4,... are unstable relative to phase 2. 

• If 63 = -1, 64 = 4, and -6 < 6 5 < -4, then 

t~\h-f 2 ) ~ (6 5 + 4)t 5 <0, 
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(2.58) t- 2n (f n+1 - f n ) ~ (& 5 + 6)t 5 > , n > 3 , 

and phase 3 is stable relative to phase 2 and to phases 4, 5, ... . 

• If 63 = -1, 64 = 4, 65 = -6, and y < h < y, then 

t- 4 (/ 3 -/ 2 )^-2t 5 <0, 

r 6 (/ 4 -/ 3 )^(& 6 -§)* 6 <0, 

(2.59) r 2 "(/n+i - /„) ^ (6e - f> 6 > , n > 4 , 

and phase 4 is stable relative to phases 2, 3 and to phases 5,6,.... 

• If 63 = —1, 64 = 4, 65 = —6, 6 6 = y, and —53 < b 7 < —51, then 

t- 2n (f n+ i - fn) ^ -2t n+3 < , 2 < n < 3 , 
r 8 (/ 5 -/ 4 ) ^ (6 7 + 51)t 7 < 0, 

(2.60) r 2n (f n+l - f n ) ~ (6 7 + 53)t 7 > , n > 5 , 

and phase 5 is stable relative to phases 2, 3, 4 and to phases 6, 7, ... . 

• If 63 = -1, 64 = 4, 65 = -6, 6 6 = f , b 7 = -53, and 160 < b 8 < 162, 
then 

t- 2n (f n+ i ~ fn) ^ -2r +3 < , 2 < n < 4 , 
^ 10 (/ 6 -/ 5 ) ^ (6 8 -162)t 8 <0, 

(2.61) t- 2n (f n+l - f n ) ~ (6 8 - 160)t 8 > , n > 6 , 

and phase 6 is stable relative to phases 2, 3, 4, 5 and to phases 7, 8, 

• There exists B 9 as follows. If 63 = —1, 64 = 4, 65 = —6, b§ = y, 6 7 = 
-53, 6 8 = 160 and B 9 < b 9 < B 9 + 2, then 

t- 2n (f n+ i - /„) ^ -2t n+3 < , 2 < n < 5 , 
r 10 (/ 7 -/ 6 ) ~ (b 9 -B 9 -2)t 9 <0, 

(2.62) t- 2n (f n+1 - f n ) ~ (6g - 5 9 )t 9 > , n > 7 , 

and phase 7 is stable relative to phases 2, 3, 4, 5, 6 and to phases 8, 9, . . . 
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3. Phases 0, 1, 2 
For n — 0, (12. 3ft takes the form 

(3.1) ^ = ^'£11^(7) 

{7} 7 

with 

(3.2) V(t) =t5l7l-l7n/n| e ^| 7 n{ 2 =i}| 



so that 



'■ T (o;l 



log(Z A ) = M |A 1 | + £V 3 

a; 

(3.3) =«|A 1 |+£V T H+ £ 



while 



iog(^) = E ^» + E + E 
(3.4) = e + E (") + E 



Therefore 
(3.5) 

=u|Aii+ E V (<*>)- E v t h+ E v#m- E ^i( w ). 

giving 

/ t - /„ = u + ( e ~ u t 2 + 2e- 2u t 3 ) - {t 2 + e"t 2 + 2t 3 + 2e 2u t 3 ) + t 3 + 0{t A ) 

(3.6) = (63 - l)t 3 + 0(t 4 ) 

For n = 1, in order to use (I2.12p . we need Ai(u) , A 2 (u) , B^u) , B 2 (0) ■ The 
expansion 

(3.7) fA^u) = e u t 2 + 2e 2u t 3 + 6e 3 V + e 4u t 4 - e u t 4 - |e 2u t 4 - 2e 2 V + C(t 5 ) 
gives 

£ 2 (Ai(u) - Ai(0)) = (e u - l)t 2 + 2(e 2u - l)t 3 + 6(e 3u - l)t 4 + (e 4u - l)t 4 

- (e u - l)t 4 - \{e 2u - l)t 4 - 2(e 2u - l)t 4 + 0(t 7 ) 

(3.8) = t 4 + (63 + 4)t 5 + (64 + 46 3 + 17)t 6 + 0(t 7 ). 
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We then compute A 2 {u) using P 2 , Q 2 : 

P 2 = (1 + 4t - t - At 2 ) + (12t 2 + 6t 2 + At 2 - 6t 2 ) + 0{t 3 ) 

(3.9) = 1 + 3t + 12t 2 + 0(t 3 ) 

where the first parenthesis is adapted from Fig. 2 and the second from Fig. 3. 
Also, adapted from Fig. 7, 

Q 2 = 4t 2 - At 3 + 12t 3 - 5t 3 + 2t 2 + 12t 3 + 0{t A ) 

(3.10) = 6t 2 + 15t 3 + C(t 4 ) 

giving 

t A A 2 (u) = e¥(l + 3t + 12t 2 ) + Qe 2u t 6 + 0(t 7 ) 

(3.11) = t 4 + 3t 5 + 19t 6 + 0{t 7 ) 

and 

(3.12) fA^u) - t 2 Ai(0) - t 4 A 2 (u) = (b 3 + l)t 5 + (6 4 + Ab 3 - 2)t 6 + 0(t 7 ). 
Then 

(3.13) t 2 5i(w) = e u t 3 + 2e 2u t 5 - e u t 5 + 0{t 7 ), 

(3.14) t A B 2 (0) = t 5 + 0(t 7 ), 

(3.15) t 2 B 1 (u) - t 2 B 1 (0) - t 4 B 2 {0) = b 3 t 6 + 0(t 7 ) 
so that finally 

(3.16) f 2 - h = (b 3 + l)t 5 + (6 4 + 36 3 - 2)t 6 + 0(t 7 ), 
which completes the derivation of (jl.6p . 
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4. Recursion diagrams, n > 3 

For the recursion relations (12.17ft relating n to n + 1, we consider ways in 
which a cluster uj E I n ,W can be extended to produce a new a/ e I n+ i,W. 
One choice is that one or more polymers in uj may be extended without adding 
polymers or changing incompatibility relations within uj. Then the combina- 
toric factor in (12.61) is unchanged, only the ^(7) for the extended polymers 
change, and it remains to find a geometric factor, the number of ways to 
extend the polymer, or the number of diagrams of a given type. 

Next, one may have uj' = uj U {7'} with the new polymer incompatible with 
only one polymer from uj. Then ( 12. 6ft gives ip T {cu') = — Lp(^f')(p T (uj), with <p T (uj) 
taking into account possible polymer extensions as in the first case. 

Next, one may have uj' = ujU {7'} with the new polymer incompatible with 
two polymers 71,72 from uj. At the order considered here, one may assume 
that 71 7^ 72 and that uj = {71,72} or uj = {70,71,72}- Then f |2.6j) gives 



(4.1) = -2^(7 , 



with <£> t (uj) taking into account possible polymer extensions as in the first case. 
Formula (14.11) occurs in the 2nd and 3rd diagrams in the 2nd line for 
and in the 3rd and 5th diagrams in the 2nd line for Qn+i- 




-6t 3 Q 2 n a -2t 3 Q 2 n a -8t 3 Q 2 n b +2tR^ +8t 2 R n +t 2 R n +tR 2 n + ... 

Fig. 10. Recursion for Q\ +1 - 
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Next, one may have u' = uj U {7^, 7 2 } with each of 7^, 7 2 incompatible with 
at most one polymer in w. If 7 2 7^ j[ to and 7 2 ~ lu, or 7 2 9^ u; 7^ 7J and 
7 2 ~ t(, then (12.61) gives (f T (u)') = (p('y f 1 )(p('y' 2 )tp T (uj), with <£> T (w) taking into 
account possible polymer extensions as in the first case. If 7 2 7^ j[ 7^ uj 7^ 7 2 
and 7 2 7^ 7i, and 7^ and 7 2 are incompatible with the same polymer in uj, then 
dH gives 

(4.2) if (to') = 2 ( p(j l )<p(j l ) ( p r (u) , 

with ip T (uj) taking into account possible polymer extensions as in the first 
case. Formula (14. 2 j) occurs in the 5th diagram in the 1st line for Q^+i and in 
the 5th diagram for R^ +1 and in the last diagram for -R^ +1 . If 7 2 = j[, then 
ip T (u)') = ip(j' 1 )(p('j2)(p T (uj), with tp T {uj) taking into account possible polymer 
extensions as in the first case. 
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Rn+l — 18i 4 P n 



l8t 5 P n 



^Q 1 



-2At 4 Q % 



-t 2 R n + . 




Rn+l — 32t^ P n 



-m 5 p n 



-28i 6 P„ 



-12t b P„ 



+8t b P n 



-8t*Q l n 




R n +i — 12i P n 



+6t b P n 



+5t b P n 



+8t 6 P n 



+ . 



S n+1 = 4t 5 P n 



-m 6 P n 



+ 16t 6 P n 



+28t 6 P n 



+4i 6 p„ 



-2t 4 Q 7 



Fig. 12. Recursions for R\ +1 , R-n+i, P„+i an d S n +i- 



Factors larger than the ±2 in (14. ip and (14. 2ft are possible for extensions 
uj' . At the given orders, though, such factors do not appear in our formulas 
for ip T {oj') or contribute to the recursion formulas (I2.17p . because the added 
polymers, 7', or j[ and 7 2 , do not create new cycles in the incompatibility 
graph other than possibly cycles of length 3, namely 7' 7^ 71 7^ 72 7^ 7' or 
il / 7 t 7 2 / 
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5. Diagrams for the 7/8 transition line 















(i) 






















c 


A 

J 


) -4("a 

3 




+1 


KV)0 


% — 


2) +22("; 

□ (2) 


i 


+4(«" 


2 ) 


+64 (V) 

r - 


- 1 

r - - - 


6(n 



-64 (V) +10 +10 (V) +4(n-2) -16 (V) ^(Y) 

+8(n-l) -2(n-2) -2(n - 2) 




+176n +32n +8n -124n -32n +32n 




□ 



+54n -20n -20n -An -48n +40n +20n 

Fig. 13. P n : i 4 terms dependent on n. Continuation down- 
ward from levels containing two cubes, as in configuration 
(1), may be from below either cube. Configurations (2) 
are excluded from the preceding two diagrams. For (3) see 
Fig. 14. 
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X X 
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Fig. 14. Terms contributing to (3) in Fig. 13. Top view; 
x represents a possible location of the column below. 



□ 



(V) +8(V)(™- 2 ) +4( n 2 1 ) +64 (V) -16(n-3) +32(n - 2) + 




-16 (V) -4(n-2) 



+!0(V) 
-2(n-2) 



+10 +4(n-2) 
-2(n-2) 



□ 



+8(n-4) +5(n-l)-l +1 -("2 1 ) -8(n - 2) 

Fig. 15. P„: i 5 terms dependent on n. 
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